Bs Semileptonic Decays to Ds and D* in Bethe-Salpeter Method* 
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Using the relativistic Bethe-Salpeter method, the electron energy spectrum and the semileptonic 
decay widths of B° -> D~l+vi and B° Dl'l+ui are calculated. We obtained large branching 
ratios, Br{Bs DsCVe) = (2.85 ± 0.35)% and Br{Bs -5> DleVe) = (7.09 ± 0.88)%, which can be 
easily detected in the future experiment. 
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The light i?" and meson decays have been precisely 
measured in experiment and well studied in theories 
Unlike and B'^ mesons, the heavier Bs meson can- 
not be studied at the B-factories operating at the T(4S') 
resonance. It can be produced at the T(5S') resonance, 
Belle collaboration has collected large number of Bs sam- 
ples in the year 2006 and 2008, and they have presented a 
preliminary result of the inclusive semileptonic branching 
fraction 0: Br{B° X+e'v) = (10.9±1.0±0.9)%, and 
in the forthcoming LHCb experiment, besides 5° and B'^ 
mesons, an abundant number of Bs meson will be pro- 
duced, so we have been provided a chance to study the 
properties of Bs meson and its various decay channels. 

The physics of Bs has become hot topic in recent years, 
for example, the non-lcptonic Bgdecays have been exten- 
sively studied in the literature [3-01 • Whereas the inves- 
tigation of semileptonic decays of Bs to heavy meson is 
relatively modest, for example, though there are the mea- 
surement of Bs — > DsXl'^ve [§], the pure decay modes 
of Bs — Dsl^vi and Bs — >■ D^i'^vi have not been de- 
tected separately in experiment. The semileptonic weak 
decays of heavy flavored mesons are certainly interest- 
ing, it not only yield some of the most useful informa- 
tion on the elements of the CKM matrix, CP-violation 
and flavor-violation, but also increase our opportunities 
of discovering new physics. 

In this letter, we will study the exclusive semilep- 
tonic decays of Bs to Dg and 2?* mesons based on the 
instantaneous relativistic Bethe-Salpeter (BS) equation 
(Salpeter equation). Since the BS equation is a full rel- 
ativistic equation to describe bound state, so we choose 
this method to consider the relativistic corrections in Bs 
decays. Another important thing in this decay is that the 
Bs is much heavier than the final state, there are large 
recoil effect, this should be dealt with correctly. In our 
calculation, besides the relativistic dynamic BS equation, 
we use Mandelstam formalism to calculate the tran- 
sition amplitude with invariant variables as input, so we 
can correctly consider the recoil effect, and give reliable 
calculations. 
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Fig. 1. Feynman diagram corresponding to the 
semileptonic decays — > X -|- -|- 

For the semileptonic decays Bs X + ^ vg (here 
X denote as Ds and I?*) shown in Fig.l, the T-matrix 
element is: 



^^)v,{X(P',t)\J^\Bs(P)) (1) 



where Vij is the CKM matrix element, is the charged 
current responsible for the decays, P, P' are the momenta 
of the initial state Bs and the final state X, respectively. 

To evaluate the exclusive semileptonic differential de- 
cay rates of Bg meson, one needs to calculate the 
hadron matrix element of the weak current sand- 
wiched by the Bs meson state as the initial state and 
a single- hadron state of the concerned final state, i.e., 
(X(P',e)|Jp|Ps(P)). It is well known that the Mandel- 
stam formalism 10] is one of proper approaches to com- 
pute the hadron matrix elements sandwiched by the BS 
wave functions of the two bound-state, no matter how 
great the recoil momentum carried by the elements will 
be. With the help of this method and the instantaneous 
approximation [ll|, the hadron matrix elements in the 
center of mass system of initial meson can be written as 
[TMl: 



{X{P')\J^\Bs{P)) 
dq 



(27r) 



rTr 



(2) 



where M is the mass of the initial meson Bg , q and q' are 
the relative three-momentum of the quark-anti-quark in 

-P',P' 
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the meson Bs and X, respectively, q' = q > , , 

is the three dimensional momentum of finial hadron state 
X , is the component of BS wave function projected 
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onto the positive energy for the relevant mesons, and 

Considering the semileptonic decay Bg ^ D^+i^ + vi 
as an example, for the initial state pseudoscalar meson 
(J^ = 0~) Bs, the positive energy wave function takes 
the general form [Tsj : 
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(3) 



where — (0,q) and = ^777,f + q^, fi{q) are eigen- 
value wave functions which can be obtained by solving 
the full Salpeter equations. As for the finial state vector 
meson (J^ = 1 ) D*, the positive energy wave function 
takes the general form 16|: 



1 r 



q'^-e'^{E + F + G r' 4'J] (4) 

where e is the polarization vector of meson, and 
A, B, C, D, E, F, G, H are defined as: 



A = M' 
B 
C 
D 



h{q')-h{q') 



7772 W J + m'iU)'2 



Mq')~h{q') 



uj'iLj'2 + m[jn2 + q' 

h[q')^k{q')^ 
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M'{uj[uj'2 +m[m'2 - q' ) 

7n\ + 777o 



M'^ h{q')-h{q')- 



F 



h{q) + h{q') ]^ ^ 

LOl + UJ2 



M'^{uj[uj!2 + m[m'2 - q'' 
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Mq')-Mq') 



h{q') + h{q')- 



G = ^ 



1 

M' 



h{q') + h{q') 



V6{q')M' 
m2Uj'i + 777']^ 



H = 



1 

JjJ2 
2/5(9') 



'hiq')4i^ + hiq') 



(5) 



(777'j^ + 7)7,2 )(a;5.a;2 + Tn'^m'2 + q'^) 
where M' is the mass of the final meson Z?* and E' = 

Now, using Eq. (2) and the formula as follows @: 
1 1 



8A-/ (27r) 



iTl^dEfdE' 



(6) 



the concerned widths of the semileptonic decays are cal- 
culated out finally. 
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Fig. 2. The electron energy spectrum corresponding to the 
semileptonic decays Bs — > DsS'^Ue 
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Fig. 3. The electron energy spectrum corresponding to the 
semileptonic decays Bs D'e'^Ue 

In the numerical calculation, there are some parameters 
have to be fixed. The input parameters for the masses of 
quarks are chosen as follows: a — e = 2.7183, a = 0.06 
GeV, A = 0.1 GeV2, Aqcd = 0.15 GeV and 777,, = 4.96 
GeV, 7?7c = 1.62 GeV, nis = 0.5 GeV. For pseudoscalar 
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mesons Bg, Dg and vector meson D*, we choose Vo = 
-0.194 GeV, Vo = -0.34 GeV and Vq = -0.15 GeV. 
With the parameters, we obtained the masses: Mb^ = 
5.3663 GeV, Md, ^ 1.9685 GeV and Md- = 2.1123 
GeV. The value of the CKM matrix elements we used in 
this paper is: Vcb = 0.0412 [§]. 

In Fig. 2 and Fig. 3 we show the spectra of the electron 
energy for the decays Bs — DsC^Vf. and Bs — > D*e^i^e, 
respectively. 

Our numerical results of the exclusive semileptonic de- 
cay widths for Dg and D* final state are: 

r(B, ^ Dgci^e) = (1.27 ±0.15) X 10"" GeV 
r{Bg Dleve) = (3.17 ±0.39) x 10"" GeV (7) 

Correspondingly, the branching ration of the two de- 
cays are: 

Br{Bg Dscve) = (2.85 ± 0.35)% 

Br{Bg Dlei^e) = (7.09 ± 0.88)% (8) 

so these two channels have large branching ratios, and 
with two electronic particles in the final states, these two 
modes are easy detected in experiment. Our results are 
close to the existing theoretical results by Blasi et al [13] , 
which are r{Bs DsCVe) = (1.35 ± 0.21) x 10"" GeV 
and T{Bg D^ei^e) = (2.5 ± 0.1) x 10"" GeV. 



We also use this method to calculate the Bg semilep- 
tonic decay to B and B*, though these modes have a 
favor CKM matrix element, Vus = 0.2255 0, they are 
deeply suppressed by phase space. And the decay widths 
are: 

r{Bs Bei^e) = 0.85 x 10"^" GeV 

r{Bg B*eve) = 1-04 X 10"^^ GeV (9) 

with the corresponding branching ratios: 

Br{Bs Beve) = 0.19 x 10"^ 

Br{Bs B*eve) = 0.23 x 10"** (10) 

which are too small to be measured. 

In conclusion, we have calculated the decay widths 
of the exclusive semileptonic Bg decays to and D* 
mesons by means of the instantaneous BS equation 
method. We find these two modes have large branching 
ratios up to 10% of the full decay width, which consist 
with the preliminary result of Belle. But the CKM ma- 
trix element favored channels, B^ semileptonic decays to 
B and B* ^ are suppressed deeply by the phase space and 
can not be measured easily. 
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